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Abstract 

Zero-field specific heat of the single crystals i?Ni2B2C {R— Er, Ho, Dy, Tb, Gd) was measured within the 
temperature range 0.1 K< T <25 K. Linearized spin wave analysis was successfully applied to account for 
and to rationalize the thermal evolution of the low-temperature magnetic specific heats of all the studied 
compounds (as well as the one reported for TmNi2B2C) in terms of only two parameters, namely an energy 
gap A and a characteristic temperature 0. The evolution of and A across the studied compounds correlates 
very well with the known magnetic properties. 0, as a measure of the effective RKKY exchange couplings, 
scales reasonably well with the de Gennes factor. A, on the other hand, reflects predominately the anisotropic 
properties: - 2 K for GdNizBzC, ~ 6 K for ErNiaBaC, - 7 K for TbNizBzC, and - 8 K for each of 
HoNi2B2C and DyNi2B2C. The equality in A of HoNiaBaC and DyNi2B2C, coupled with the similarity in 
their magnetic configurations, indicates that a variation of x in the solid solution Ho2,Dyi_2;Ni2B2C {x < 0.8 
and Tc < T/v) would not lead to any softening of A. This supports the hypothesis of Cho et al (PRL 
77,163(1996)) concerning the influence of the collective magnetic excitations on the superconducting state. 
This work underlines the importance of spin-wave excitations for a valid description of low-temperature 
thermodynamics of borocarbides. 

PACS numbers: 

I. INTRODUCTION 

Wide varieties of magnetic ground structures are manifested in i?Ni2B2C series (see e.g. Refs^*^ 
and Tab. I): members with R =Ho, Dy, Nd, Pr adopt commensurate antiferromagnetic (AF) ar- 
rangements while those with R =Tm, Er, Tb, Gd, assume squared-up spin-density wave (SDW) 
states. These ground structures are stabilized by a fine balance of exchange, crystalline electric field 
(CEF), and magnetoelastic forces. Under magnetic field or temperature variation, most of these 
structures undergo a cascade of phase transformations, yielding rich varieties of field-temperature 
(H — T) phase diagrams ( see e.g. Refsf^*^*^). Interestingly, most of the zero-field magnetic 
states of R=Tm, Er, Ho, Dy coexist with superconductivity, presenting model compounds wherein 
the interplay between superconductivity and magnetism can be investigated. Such investigations 
revealed that the superconductivity, though much influenced by, has a very weak influence on the 
prevailing magnetic order:^ the energy gain due to the onset of magnetic order dominates by two 
order of magnitudes over that due to the onset of superconductivity. 

The zero-field part of the H — T phase diagram of HoNi2B2C is particularly interesting: su- 
perconductivity sets-in at Tc ~ 8 K. Just below this point, an incommensurate spiral state, 
kc ~.92c, develops (see e.g. Refs4*i£). Furthermore, at ~6.3 K, an additional modulated state 
with ka —0.55a emerges and around 5 K a deep minimum in Hc2 develops.^ At T^r =5 K, an 
orthorhombic lattice distortion^'^ sets-in and, concomitantly, both the spiral and the a-axis modu- 
lated states are replaced by a commensurate AF structure which coexists with superconductivity 
down to the lowest measurable temperatures. Generalized susceptibility calculatioi>ii related these 
modulated state to maxima in the exchange coupling transform, J{k). 

Remarkably, the magnetic ground structures of the heavy and magnetic i?Ni2B2C compounds 
(see Tab. I) are particularly simple: in spite of the manifestation of an orthorhombic distortion and 
a liquid-helium-temperature magnetic modulation, their ground structures are either an equal- 
amplitude, AF-type squared-up state (as in R =Tm, Er, Tb, Gd) or an equal-amplitude, colinear, 
commensurate AF state (as in R =Ho, Dy). Then, it is of interest to investigate whether the 
low-temperature thermodynamics of these ii!Ni2B2C can be described in terms of small-amplitude 
spin-wave excitations and, in addition, to elucidate the character and dimensionality of these 
excitations. These magnetic excitations can be probed by various techniques, among which is the 
magnetic specific heat. We carried out extensive zero-field specific heat measurements on five single 
crystals i?Ni2B2C (i?= Er, Ho, Dy, Tb, Gd) covering at least the range 0.1 K< T < Tat. These 
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specific heats, together with that of TmNi2B2C (Refj^^i), reveal a diversified and wide varieties 
of thermal evolutions. None the less, based on a phenomenological model, all of the specific heat 
curves can be systematized in terms of only two parameters, namely an effective exchange coupling 
and a magnetic anisotropy interaction. 

The format of this paper is as follows: In Sec. II, we derive an approximate, but of wide ap- 
plicability, expression for the magnon specific heat. Experimental techniques and procedures are 
described in Sec. III. Results and their analysis are described in Sec. IV and discussed in Sec.V. 



II. APPROXIMATE EXPRESSIONS FOR MAGNON SPECIFIC HEAT OF i?Nl2B2C. 



The magnetic structures of i2Ni2B2C (i?=Tm-Gd) can be visualized as magnetic layers that 
are stacked along the c-axisii The most dominant interactions are the RKKY and anisotropic 
couplings. The former can be approximated by effective isotropic couplings while the latter (a 
combination of dominant CEF and weaker dipolar and anisotropic exchange forces) by an easy- 
axis anisotropy field this Ha representation is convenient for spin- wave calculation and is 
valid for the low-temperature orthorhombic-distorted equal-amplitude collinear AF or squared-up 
phases of borocarbides (see Sec.V). Considering the above mentioned magnetic arrangement and 
the crystal symmetry, the magnetic couplings can be conveniently divided into two classes: 
that couples moments i and j within the same layer (denoted as A or B) and J^^ that couples 
moments from different layers. Then at zero external field, the following Hamiltonian is expected 
to capture most of their low-temperature properties: 

H = - ^ j[jSi.Sj+ ^ J^^Sf.Sf - gfisHa'^Sf + gf-isHa-'^Sf (1) 

<ij>,L£A,B <ij>A,B jeA j&B 

The first and second sum represent, respectively, interactions within the same layer and among 
different layers. The last two terms represent the anisotropic interactions. By standard spin-wave 
analysis, we obtained the following dispersion relation: 



nwk = yj [SJi. (0) - SJ^ {k) + 5 J|| (0) + g^XBHa] ' - [5 J|| {k)] ' (2) 

where J{k) = ^ Jij ex.p{ik{ri — Vj)). ( J_|_) represents the Fourier-transform parallel (perpendic- 
ular) to the c-axis. The energy gap, A = fkJk=Oj is: 



A = ^JiofiBHaf + 2J\\{Q)gfiBSHa (3) 

which defines an AF resonance frequency similar to the uniform mode of ordinary AFs(i^*i^'i^*i^ 
Evidently, (i) A is zero whenever there is no anisotropy and (ii) A does not depend on the type 
nor the strength of the intralayer coupling. 

For evaluating the magnon specific heat (or other thermodynamical quantities), an explicit 
expression of J{k) is required. In the absence of such an expression and for low-temerpature 
range, it is a common practice to assume a long-wave limit. Here, we restricted the expansion of 
J{k) to the nearest neighbors only, leading to: 



huk « ^A2 + C^kl + Cykl + C,kl (4) 

where Cx {~ Cy for weak orthorhombic distortion) and Cz are functions of the exchange couplings 
and geometrical factors (a and c are unit-cell parameters): 

Cx = 16(Jx + J\\)J\\S'^a'^ + 2J±S{gfiBHa)a'^ (5) 
a = 16 J^S^c^ 



3 



Then, the zero-field magnon specific heat is (rewritten so as to conform with the notation of Refj^): 

oo 

Cm{T) = 3^R {A/ef {A/2Ttt^) ^ [BessdK{2, mA/T) + BesselK{4, mA/T)] (6) 

m=l 

where BesselK represents the modified Bessel function and 

e = z\J,ff\S = ll 2,i2{C^CyC-;)-2 /{a^c) (7) 

is a characteristic temperature, based on which | Je//| can be defined as being an effective exchange 
interaction that couples the magnetic moment to its z nearest neighbors. Notice that both EqEl 
and EqI7| account as well for the weakly orthrohmic-distorted state. 
For T < A, EqEl reduces to the exponential form: 

3 „ . 7 

CM(r)=.2^— ^exp(-A/r) (8) 
while for isotropic compounds or T ^ A, it reduces to the high-temperature limit>i^*i^ 

CM{T)^^-^R{T/ef (9) 

in this case, Eqd highlight the useful definition of \ Jeff\- 

A long-wave dispersion relation for an isotropic quasi-Sd case can be derived from Eq|ll if we 
set Cz^Cx (| J|| I <C J±) and Ha = 0: 

hwk^SJ\\S + J^Sa^ikl + kl) (10) 

Then, to lower order in SJ\\S/T, one obtains Cm{T) = ttRT /12SJ±_ which reproduces the leading 
linear-in-T term in the expression reported by Movshovich et al who (starting from a quadratic 
dispersion relation and including correction for the 2d and the magnon-magnon interaction) ob- 
tained for the range | Jy | < T < Tjv: 

Cm{T) = {7:R/12){T/SJ± - QJ\\/tt'^J^ + 4J|| S/Svr^ J^T) (11) 

It is worth remarking that Eqs lll4l6l are of a more wide applicability than our above analysis 
might have suggested. Furthermore, depending on only 9 and A, a variety of expressions for the 
magnetic specific heat can be derived; each admitting different limits for its thermal evolution 
depending on the relations among T and A (compare Eql^Jwith the limit Eas l8l9irT|l . Based on 
such a scheme, one is capable of rationalizing the vast varieties of the low-temperature thermal 
evolution of thermodynamical quantities (such as Cm{T)) encountered in these (and any series 
similar to) borocarbides. It is reminded that this analysis is not adequate for the description of 
the contribution of the modulated state since the involved ordered components do not have equal 
magnitudes. Finally, the above zero- field treatment is obviously not appropriate for the description 
of the field- induced metamagnetic phases. 



III. EXPERIMENTAL 

Single crystals of i?Ni2B2C (i?=Er-Gd) were grown by floating zone method»ii Physical char- 
acterization were carried out utilizing a dc magnetometer [//<50 kOe, 2 K <r<20 K], a mutual- 
induction ac-susceptometer [z^=250-500 Hz, Hpp ssl Oe, 30 mK <T<20 K], and a dc resistivity 
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set-up [H<50 kOe, 2 K <T<20 K]. Structural, magnetic and transport characterizations are in 
agreement with pubHshed results. 

The temperature-dependent specific heat was measured on a semi-adiabatic calorimeter [80 mK 
<r<25 K, precision better than 4%]. The total specific heat curves measured above 2 K are in 
agreement with the reported dataiS*iSii2i2ii However, we observed some discrepancy between the 
absolute values of Cm{T) of single crystal and polycrystalline samples: though both specific heats 
were found to be given by approximately the same functional form, the absolute values of the fit 
parameters [9 and A) differ by as much as 40%. 

For each compound, the total specific heat, Ctot, was analyzed as a sum of an electronic Ce {Cg 
when superconductivity is to be emphasized), a Debye Cd (= PT"^), a nuclear Cn, and a magnetic 
contribution Cm from the only magnetically active i?-sublattice. At temperatures of interest, Ce 
and Cd were estimated based on our specific heat characterization^^ of single crystal YNi2B2C 
(7 = 17.5 mJ/moleK^ and /? = 0.12 mJ/moleK^) which had been synthesized by the very same 
procedures as the one used for the other single crystals. Thus for each compound, 7 was taken to 
be the same value as that of YNi2B2C while /? was evaluated by normalizing the /? of YNi2B2C 
using the approximation /3 oc ma (see Tab. II). 

Within the superconducting region, C5 was evaluated as 3jT^/T^. At any rate, for all the 
studied compounds, Cm{T) is much larger than the sum of Ce and Cd- Consequently, even if 
Ce and Cd are taken as the bare values of YNi2B2C, Cm{T) would not be noticeably modified, 
ensuring that our conclusions would not be influenced. 

Cm{T) of -RNi2B2C, when available, is of dominant importance only at very low temperatures 
and was evaluated by least-square fit using the appropriate hyperfine Hamiltonian; the obtained 
parameters (shown in Tab. II) compare favorably with those of the corresponding R-meta^^ and 
i?Co2B2C isomorphsi^ 

IV. RESULTS AND ANALYSIS 

From the general feature of Cm{T) curves (Figs. 1-3, 5-6), one distinguishes four temperature 
regions: (i) a paramagnetic region, T > T/v, wherein Cm{T) is predominantly due to a change in 
the population of the CEF levels, (ii) a critical region, T ^ T/v, wherein Cm{T) is related to the 
critical phenomena, (iii) an intermediate region wherein Cm{T) reflects the magnetic character of 
the spiral/modulated states, and (iv) the low-temperature AF/squared-up states, of prime interest 
to this work, wherein the measured Cm(T) is to be confronted with Eqs Killll and therefrom 6 and 
A are to be extracted. 

Before we discuss the features of Cm{T) for each compound, a word of caution is in order: 
just as in the case of i?-metals(^ the propagation of errors due to successive subtraction of Ce(T), 
Cd{T) and Cn{T) would eventually influence the absolute value of Cm{T). 

A. GdNiaBaC 

Below T/v, the zero- field magnetic structure^ is a transverse sine- modulated type with ka that 
changes from 0.551a at T/v to 0.550a at T/j, where a spin reorientation process sets-in. Below Tr, 

— ^ * 

ka reverts course and increases monotonically till reaching 0.553a at 3.5 K. 

Cm(T) of GdNi2B2C (Fig.l) reveals the onset of the magnetic order at T/v =19.5 K and the spin 
reorientation process at Tr =13.5 K, in agreement with earlier studies i^'^S'^SiiLS^ The thermal evo- 
lution of Cm(T) within the amplitude-modulated state is distinctly different from that within the 
equal-amplitude, low-temperature state (see below). As mentioned above, within the modulated 
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region the linearized spin-wave analysis is not applicable and one should resort to the findings of 
Schmitt and co-workers^S. In particular, these authors demonstrated that Cm{Tn) of such a state 
suffers a strong reduction (almost 1/3) in comparison with the value (20.15 J/moleK) expected for 
an equal-amplitude AF state (our results are in excellent agreement with this statement). 




T(K) 

FIG. 1: log-log plot of Ctot{T) (triangles), Ce(T) + Cd(T) (dotted), and Cm(T) (circles) curves of single 
crystal GdNiaBaC. The continuous line represents EqEI with 6 =12.5±0.2 K and A = 1.9±0.3 K. For 
A < T < 4 K, Cm(T) follows O.OSST^ J/moleK (dashed line) which is the high-T Hmit of Eqini(see text). 

On the other hand, Fig.l shows that for temperatures below 3.5 K, Cm{T) follows faithfully Eq|Sl 
with ^=12.5ib0.2K and A=1.9ib0.3 K. The numerical value of 6 (for A see Sec.V) is physically 
acceptable as can be seen from the following arguments. First, the substitution of this 9 into 
Eq|3 yields Jeff= 0.58ib0.2 K, which is in close agreement with the value reported for HoNi2B2C 
(Refpi^) and TmNi2B2C (RefF^). Secondly, the substitution of 9 into EqEl predicts correctly 
the high temperature limit, namely Cm{T > A) =0.058r'^ J/moleK (see Fig.l). Thirdly, the 
substitution of 9 into the Molecular Field relation)^^ 

T^ = l0(^S + l) (12) 

gives Tjv = 18.8±0.3 K which is in reasonable agreement with the experimentally determined value 
of T^. 

B. TbNiaBaC 

A longitudinal SDW, accompanied by an orthorhombic distortion, sets- in at T]\fM^ The mag- 
nitude of the modulation vector decreases from 0.551o near Tjv to 0.545O at 2.3 Ki^ A weak 
ferromagnetic component develops below T]yF ^ 8 K and at lower temperature a squaring-up of 
the modulated state was observed 
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Cm{T) of single crystal TbNi2B2C (Fig.2) shows the magnetic ordering at T/v = 14.5 K and 
the WF-associated anomaly that peaks around 5.5 K. These features are in agreement with those 
reported by Tommy et aL— Similar to the case of GdNi2B2C, no attempt was made to analyze 
Cm{T) within the amplitude modulated state spanning the range 5 K< T < Tjy. Below 5 K, 
where the orthorhombic-distorted squared-up state is expected, Cm{T) follows convincingly the 
prediction of EqElwith 6^=21. 5±0. 2 K and A=7.0±0.5 K. 



FIG. 2: log-log plot of the Ctot{T) (triangle), Ce{T) + CD{T) (dotted), Cn[T) (dashed) and Cm{T) (circle) 
of single crystal TbNi2B2C. The continuous line represents EqElwith A = 7.0±0.5 K and = 21.5±0.2 K. 



C. DyNizBzC 

This compound develops a commensurate AF structure below T^r with moments arranged in 
an identical manner as that of HoNi2B2Cfi*^ Moreover, superconductivity coexists with this AF 
order below Tc ~6 K. In contrast to other Ni-based AF superconducting borocarbides, DyNi2B2C 
presents the following distinct features: (i) no zero- field incommensurate or modulated state was 
reportedji*^ However for T < 2 K, anomalously large hysteresis and pronounced reentrant effects 
were observed for the field range 1 kOe< H <5.3 kOe(24 (ii) the superconductivity emerges within 
a well developed AF order (Tc < T^r) and that Tc is extremely sensitive to nonmagnetic doping. 

Ctot{T) (see Fig. 3) reveals the onset of the AF order at T/v = 9.5 ± 0.2 K. Within the accuracy 
of our measurement, the superconducting jump at ~6 K is too small to be resolved. On carrying 
out the analysis of Ctot into its components (Cg, Cd-,Cn-, and Cm), we observed an anomalous 
contribution peaking at 1.2 K and having features reminiscent of a Schottky-like contribution. 
Accordingly, it was approximated by the standard two-level relation: 




1 



10 



T(K) 



Csch{T) = R{^f exp(|)/ [l + exp(|) 



(13) 
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where 5 is the energy separation. It was found out (see the inset of Fig. 3) that 5 = 2.9 K and that 
only 0.062 molar fraction was involved. Moreover, the fit was satisfactory for the high temperature 
tail but not so good at the lower temperature part, suggesting that a multi-level Schottky contribu- 
tion might be more appropriate. However, for the present discussion, the above two-level approxi- 
mation is sufficient. It is highly possible that such a contribution is due to 6% defect /impurity which 
is beyond the limit of detection of our X-ray structural characterization. Coincidently, anomalous 
hystereses effects were observed in the magnetostriction curves that were measured within the same 
temperature range. At any case, for T > 5, both Csch{T) and Cn{T) are smaller than Cm{T) (see 
inset of Fig.3). Nevertheless, we considered Cm{T) = Ctot{T) - Cs{T) - Cd{T) - Cn{T) - Csch{T) 
(see Fig.3). 

The thermal evolution of Cm{T) is shown in a log-log plot in Fig.3 and as \n{CM{T)) versus 
1/T plot in Fig. 4. In both figures, the comparison with Eq|H| (solid line) is also presented. Evidently 
over a wide range of temperatures, Cm{T < Tc) follows excellently EqElwith 6'=19.3ib0.2K and 
A=8.3±0.3 K. 




FIG. 3: A log-log plot of C(ot(T) (triangle), Ce (T) -I- Cd (T) (dotted), Cn{T) (dashed), and Cm{T) (circles) 
of single crystal DyNi2B2C. The solid line represents EqEl (see text). The inset shows the individual 
contribution of Ctot{T) (symbol), Cn{T) (dashed), the magnetic fit (thin soUd line), Csch{T) (dash-dot, see 
Ea ll3ll . and the thick solid line is the sum of all contribution. 



D. HoNiaBaC 

Ctot{T) of single crystal IIoNi2B2C (shown in Fig. 5) reveals a cascade of three transitions that are 
usually attributed to magnetic transformations. ^■'^^ The signature of the onset of superconductivity 
is too weak to be observable in our present measurements. On the other hand, for T <5 K, Cm{T) 
follows the description of EqElwith 6l=9.7±0.2K and A=8.3±0.3 K. 
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FIG. 4: 1ii(Cm(T)) versus l/T curve of single-crystal DyNi2B2C. The data (circles) are compared to EqEI 
(solid line) giving A =8.3±0.3 K and 9 = 19.3 ±0.2 K. The inset shows on a linear scale, Cm{T) (circle) 
together with the comparison to EqEKsee text). 




T(K) 



FIG. 5: A log-log plot of C(T) versus T of single crystal HoNiaBsC showing Ctot{T) (triangle), Ce(T) + 
Cd{T) (dotted), Cm{T) (dashed), and Cm{T) (circle) contributions. The magnetic contribution (circles) 
are compared to EqO (solid line) giving A =8.3±0.3 K and 9 — 9.7±0.2 K(see text) 
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E. ErNisBzC 



Two intriguing features of the H — T phase diagram of ErNi2B2C are 111 : the onset of 
the incommensurate transversely polarized SDW state with fca=0.553o* at Tjy =5.94 K, and the 
onset of weak ferromagnetism (WF) at T\yf =2.2 K. These two events (none is able to quench 
superconductivity, Tc = 10.5 K) are well evident in Cm{T) of Fig.6. Twf-, in particular, is evident 
as a change of slope that separates two distinct thermal evolutions4S Cm{T) within the amplitude- 
modulated state Tyi/F < T < and that within the squared- up state at T < Twp- In the latter 
region, Cm{T < Twf) is weh described by EqElwith 9 =7.4±0.2 K and A =5.4±0.3 K. 




FIG. 6: A log-log plot of Ctot{T) (triangle), Ce{T) + Cd{T) (dotted), Cn{T) (dashed), and Cm{T) (circle) 
of single crystal ErNi2B2C. The solid line is a comparison to Eq. Elwith 9 =7.4±0.2 K and A = 5.4±0.3 K 
(see text). 



F. TmNiaBaC 

Superconductivity sets-in at 11 K and, below T/v = 1.52 it 0.05 K, coexists with a transversely 
polarized SDW state wherein the spins are pointing along the c-axis and the modulation vector 
is (0.093,0.093,0) piii^ At lower temperatures, the incommensurate SDW state squares upji The 
magnetic specific heat of single crystal TmNi2B2C was measured by Movshvich et ali^ and was 
shown to follow EqCUwith J± ^0.8 K, Jii «0.2 K. 



V. DISCUSSION AND CONCLUSIONS 

On a linearized spin-wave (non interacting magnon gas) approach, one is always limited to 
the low-temperature region of the ordered state which in the case of the borocarbides amounts 
to being restricted to below the liquid Helium temperatures. Within that temperature range, the 



10 



model describes very successfully the magnetic contribution of the studied compounds indicating 
that Cm can be safely associated with the gapped collective excitations that propagate within the 
orthorhombic-distorted 3d AF (-type) structure. 

The successful applicability of the model to the commensurate collinear AF ground structures 
of R =Ho, Dy is understandable. To justify its applicability to the cases of R =Tm, Er, Tb, Gd, 
it is sufficient to show that their ground structures are well squared-up and that all the moments 
do have equal amplitudes and orient (or bunch) along a specific direction. The state of ErNi2B2C 
below T]YF and that of GdNi2B2C below 3.5 K provide the best illustrations of the fulfillment of this 
requirement. The colinear, equal-amplitude, and squared-up character of the state of ErNi2B2C 
below T\Yp was elegantly revealed in the neutron diffractions studies of Choi et al*^'' and Kwano- 
Furukawa et al2^. It is remarked that the presence of weak ferromagnetism (reflected as kink that 
separate oppositely oriented domains) would hardly modify this picture since the excitation energy 
of the kink is much higher than that of the magnon. Let us now discuss the character of the state of 
GdNi2B2C below 3.5 K. A recent magnetoelastic stud5*^on single crystal GdNi2B2C demonstrated 
the presence of substantial magnetoelastic and anisotropic exchange interactions, in particular 
below T < Tr. The magnitude of the e'^ strain mode (which leads to symmetry lowering from four- 
fold tetragonal to two-fold orthorhombic) is very large and increases with decreasing temperature 
leading progressively to an orthorhombic-distorted magnetic state wherein the moments, due to 
entropy arguments, approach equal amplitudes. Based on the suggestion of Detlefs et al^, the 
magnetic structure below Tr is either a transverse modulated state with moment orientation away 
from the b-axis in the bc-plane or a modified spiral-like structure which, due to the low symmetry 
(2mm) of the [100] direction, is likely to suffer fanning or bunching (becoming more stronger as the 
temperature is decreased). Tomala et a&, on exploring these two structural possibilities, argued 
that the ^^^Gd Mossbauer spectra at 4.2K <T < Tr were better fitted with a bunched spiral-like 
state. In this Mossbauer analysis, the c-component was found to be almost equal ( ~75% at 4.2 
K) to the b-component, a result that does not reflect the 2mm-symmetry of the [100] direction. 
Considering these and the above observations together with entropy arguments, one concludes that 
the low-temperature structure is either a squared-up, equal- amplitude, and collinear state or - as 
supported by Mossbauer analysis- a strongly bunched, squared-up, and equal- amplitude state. As 
far as the magnon specific heat is concerned, the interactions in both structures can be represented 
by the Hamiltonian of EqQ (allowing for a unimportant difference in the direction of the axis of 
quantization). 

Figs. 1-6 and Tab. I demonstrated convincingly that, based on only A and 6, the diverse functional 
form of the measured Cm{T) can be systematized: when both A and 6 are large, Cm{T < A) 
reflects a magnon contribution from an anisotropic magnetic structure as in R =Er, Ho, Dy, Tb. 
When A is relatively small but 6 is large, Cm{T > A) reflects a magnetic contribution from a quasi- 
isotropic magnetic structure as in GdNi2B2C. For a weak A and J||, Cm{T) reflects a contribution 
from a quasi-2(i structure as in TmNi2B2C. 

The evolution of A and 6 across the studied compounds is reasonable. 9, on the one hand, 
reflects predominately the evolution of the de Gennes factor (see Eqs.5,7 and Tab. I) as can be 
appreciated on observing that 6 scales very well with the de Gennes factors for i?=Tm, Er, Ho, 
Dy, Tb. That the experimentally determined of GdNi2B2C is a factor of three lower than the 
one expected from de Gennes scaling may be attributed to the additional dependence of 6 on Ha 
which for GdNi2B2C is the lowest. 

A, on the other hand, reflects the combined influence of the anisotropic forces and interlayer 
coupling. This is expressed by Eq|21which for the case of, say, ErNi2B2C (considering Ha ~15 kOe 
and I J||| «0.1 K) gives a value of 4 K which is close to the experimental value. The observation 
that A is nonzero for each of the studied compounds is in agreement with the reported anisotropic 
features of the magnetic and transport properties»i^ The strong anisotropy of each of R =Er, 
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Ho, Dy, Tb is in accord with what is expected from their CEF properties. In contrast, the weak 
anisotropy observed in GdNi2B2C is most probably due to a combination of anisotropic exchange, 
dipolar, and two-ion magnetoelastic couplings. 

It is interesting to discuss one particular aspect of the interaction between magnons and su- 
perconductivity in, say, R= Ho (T^r < Tc) and Dy (Tc < Tjv): the magnon- mediated pairbreaking 
effect that manifests itself in the thermal evolution of Hc2{T < T^^Tc) [Ref42i]. Noteworthy, the 
thermal evolution of Hc2{T < T]\f,Tc) of both HoNi2B2C and DyNi2B2C [Ref>2>] are very similar 
which, considering the above-mentioned similarity in their magnetic properties, suggests that the 
involved pairbreaking effects (in particular the magnon-mediated one) are similar. This, in turn, 
suggests that the magnon characteristic (say low-energy magnon spectra) in both compounds must 
be similar. This is indeed the case: the analysis of Sec.IV.C-D showed that the energy cost for 
magnon excitation in both compounds is practically equal (A ~ 8 K see Tab. I). Therefore doping 
of Ho into DyNi2B2C (up to 80% but still Tc < Tjv) would not lead to a softening of A. Then, for 
this concentration range, there should be no variation in Hc2 and Tc even though the deGennes 
factor does vary (remember that within the antiferromagnetic state, the Abrikosov-Gorkov pair- 
breaking mechanism is not valid). This provides an additional experimental confirmation of the 
hypothesis of Cho et ai^ that the magnon spectrum of (Dyi_2:Hoa,)Ni2B2C is hardly modified for 
X <0.8. In contrast, for (Dyi_;i;Ho2.)Ni2B2C {x >0.8), the onset of superconductivity occurs within 
the paramagnetic state and consequently the Dy dopant depresses Tc linearly as expected from the 
Abrikosov-Gorkov theory. 

In summary, we were able to reveal the magnon specific heat contribution of the heavy members 
of the borocarbides and to identify the expressions that describe their thermal evolution. These 
expressions (given in terms of only two physically accepted parameters) were derived from the spin- 
wave analysis of a simple Hamiltonian that consists of effective exchange couplings and anisotropic 
interactions. We investigated as well an influence of the magnons on the superconductivity of these 
AF superconductors. 

Improvements and extension of this analysis are underway. These include, on the experimental 
side, probing the magnon contribution in single crystals of i?Ni2B2C by other (microscopic and 
macroscopic) techniques and, on the theoretical side, a better and more realistic approximation of 
J{k), CEF effects, and magnetoelastic, and anisotropic exchange forces. 
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A. Tables 



TAB. I Some zero-field parameters of selected i?Ni2B2C compomids. Superconducting Tc, 
magnetic T^y, magnetic structure, propagation wave vector, and moment direction are taken from 
Ref.i. Squaring of the modulated SDW state is taken to occur at lower temperatures. The gap, A, 
and characteristic temperature, 6, were determined from the indicated equation and figure. 6exp of 
TmNi2B2C is calculated by substituting into Eqs l5l7l the fit values (J_l=.8 K and |j|||=.2) given 
by Movshvich et al.^^. OdeG is the de Gennes scaling of 6 taking that of R =Ho as a reference. 
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TAB.II List of values of 7, /3, and 6 of flNiaBaC (i?=Gd,Tb, Er, Ho, Dy, Tm). Also given are 
the expressions used for evaluating Cs{T < Tc). The nuclear hyperfine parameters a and P (see 
Refsj"^) are also indicated: The upper (lower) values of DyNi2B2C correspond to the isotope ^^^Dy 
(^^^Dy). Data of TmNi2B2C were taken from Ref>^. 
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